Distinguishability of scalar field models of dark energy 
with time variable equation of state parameter 
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The possibility of distinguishing of scalar field models of dark energy with different Lagrangians 
and time variable equation of state parameter by available observational data is analyzed. The 
multicomponent cosmological models with the scalar field with either Klein-Gordon or Dirac-Born- 
fnfeld Lagrangians as dark energy and the monotonic decreasing and increasing equation of state 
parameters are considered. It is concluded that scalar field models of dark energy with decreasing 
and increasing EoS parameters should be distinguishable at the accuracy level of forthcoming obser- 
vational data. The Lagrangians of scalar fields could be distinguished by expected observational data 
(Planck, SDSS etc.) in the case of decreasing EoS parameter, but are practically indistinguishable 
in the case of increasing one. 
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I. INTRODUCTION 



The observational data on Supernovae of type la light 
curves, CMB anisotropy and large scale structure of the 
Universe are well fitted by the cosmological models with 
scalar field as dark energy. In numerous papers (see 
e.g. the latest reviews (Jill, books (9l4llj| and citing 
therein) the physical, cosmological and astrophysical as- 
pects of dark energy are elucidated. But the main ques- 
tion "which field from the large number of candidates is 
preferable in the light of obtained up to now observational 
data?" remains without answer. The "goodness of fit" 
of many of them, which are quite different, for the same 
dataset is so close that the problem of principal distin- 
guishability of scalar field models of dark energy arises. 
In the paper we analyze this problem for the subclass of 
minimally coupled dynamical scalar fields. 

In the fluid approach the minimal number of dark en- 
ergy parameters required for description of the dynam- 
ical and clustering properties of a scalar field is three: 
the current density parameter fide, equation of state 
(EoS) parameter w and effective sound speed cf. They 
in the general case are functions of the scale factor a 
(or redshift z = a^ 1 — 1). The last function can be 
deduced from scalar field Lagrangian when it is given. 
The function w(a) usually is defined ad hoc. At the cur- 
rent level of possibility of constraining the dark energy 
parameters by observations the simplest linear function 
w(a) = w + w a (l — a) [H, HH is widely used. Two con- 
stants wo and w a (present values of w and its first deriva- 
tive with respect to the scale factor with opposite sign) 



'Electronic address: novos@astro.franko.lviv.ua 
^Electronic address: olka@astro.franko.lviv.ua 
t Electronic address: apus@astro.franko.lviv.ua 



are determined together with density parameter fide and 
other cosmological parameters, the minimal set of which 
contains six ones: density parameter of baryons tub, den- 
sity parameter of cold dark matter uj c d m , Hubble constant 
Hq, spectral index of initial matter density power spec- 
trum n s (scalar mode), amplitude of initial matter den- 
sity power spectrum A s and reionization optical depth 
r re i . The current observational data on Supernovae type 
la light curves, CMB anisotropy and large scale structure 
of the Universe allow the possibility to determine most 
of them with high accuracy. The exception is parame- 
ter w a '- its uncertainty is so large, that the character of 
variation of EoS parameter - increasin g or decreasing - 
is not recognized (see, for example, |15l - [r7| ). 

In the paper [1 81 ] we have determined the parame- 
ters of scalar field models of dynamical dark energy 
with the other parametrization of EoS, which follows 
from assumption that the temporal derivative of dark 
energy pressure is proportional to the temporal deriva- 
tive of its energy density. This parametrization in- 
volves two free parameter wq and w e , which are the val- 
ues of EoS parameter at current and early epochs re- 
spectively. The mean likelihood distribution obtained 
by the Markov Chain Monte Carlo (MCMC) method 
for the available dataset on CMB anisotropy, power 
spectrum of spatial inhomogeneities in distribution of 
galaxies and SN la photometric curves has two peaks 
with corresponding sets of best fitting parameters = 
(fide, w , w e , uj b , uj cdm , H , n s , A s , T rei ), i = 1, 2 (Ta- 
ble [TJ . The first of them corresponds to the scalar field 
model with decreasing EoS parameter, the second one 
to the scalar field model with increasing EoS parameter. 
Since the — logL's for both p x and p 2 are close, we have 
concluded that these models are indistinguishable by the 
used dataset. This situation is the same for both classical 
and tachyonic scalar fields. Moreover, they are indistin- 
guishable too. 
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In this paper we analyze the quantitative differences of 
this four models of dark energy - classical and tachyonic 
scalar fields, each with p x and p 2 parameter sets, and 
discuss the possibility of their distinguishing in the light 
of available and expected observational data. 

The paper is organized as follows. In Section [TT] we 
discuss the scalar field models of dark energy and de- 
termination of their parameters. Section Mil is devoted 
to the analysis of background dynamics in the models 
with increasing and decreasing EoS parameters of scalar 
field and diferences of SN la light curves. In Section 
IIVI the differences of matter density and CMB tempera- 
ture fluctuations power spectra of all models are analyzed 
and compared with corresponding observational uncer- 
tainties. Section [V] summarizes the results of analyses 
and presents the conclusions. 



II. MODELS 

We analyze the dynamics of expansion and the large 
scale structure formation of the spatially flat Universe 
filled with the non-relativistic particles (cold dark matter 
and baryons), relativistic ones (thermal electromagnetic 
radiation and massless neutrino) and minimally coupled 
scalar field with given Lagrangian. We apply Einstein 
equations in the Friedmann-Robertson- Walker (FRW) 
metric for background dynamics and Einstein-Boltzmann 
system of equations in the synchronous gauge for the evo- 
lution of linear perturbations. The line-element of FRW 
metric is ds 2 = c 2 dt 2 — a 2 (t)d<j 2 , where da 2 is time- 
independent Euclidian metric of 3-space and a(t) is the 
scale factor, which we normalize to 1 at current epoch 
to: a(to) = 1. Below we put the speed of light equal to 1 
(c = 1), therefore the time has the dimension of a length, 
the Hubble constant of an inverse length, the mass has 
dimension of energy, velocity is dimensionless and so on. 
When Hubble constant appears in the traditional dimen- 
sion of [km/(s-Mpc)] the speed of light appears in [km/s]. 



A. Scalar field models 

The next specification of the scalar field model is used: 
its equation of state is Pd e (a) = w(a)pde{oi) and the tem- 
poral derivative of the dark energy pressure is propor- 
tional to the temporal derivative of its energy density, 
Pde = c 2 r^pde, where the coefficient c 2 , de ,, often called 
"adiabatic sound speed" , is constant. The last condition 
together with the differential energy-momentum conser- 
vation law, which in the Friedmann-Robertson- Walker 
metric is 

a Pde = -3(1 + w(a))p de , (1) 
gives the ordinary differential equation for w(a): 

aw' = 3(1 + w)(w-c? a{de) ), (2) 



where a prime denotes the derivative with respect to the 
scale factor a. It has the analytic solution: 



w(a) = 



(l+U7 e )(l + tOo) 



1 + w - (w - w e )a 3 ( 1+u, =) 



1, 



(3) 



where 



'a (de 



we have denoted by w e since it corresponds 



to the EoS parameter at the beginning of expansion 
(w e = w(0)), which we call early EoS parameter. One 
can see also that other constant value here, the integra- 
tion constant wq = w(l), is EoS parameter at the current 
epoch. The differential equation ([T]) with © has the an- 
alytic solution for pd e - 

Pde(a) = p de ' — — , (4) 

1 T" W e 

where p d °] = ZH^Q.^e/SirG is the dark energy density at 
current epoch. It is regular function of a for any values of 
wo and w e , excluding singular point a — for w e > — 1. 
The pressure of fields with parametrization ([3]) is regular 
function of a too: 



Pde (a) = w e p de (a) + (w Q - w e )p\ 



(0) 



(5) 



which is the generalized linear barotropic equation of 
state. That is why we call such models "the scalar field 
models of dark energy with barotropic equation of state" . 
For the quintessential scalar field (w e , wq > — 1) both 
functions, pde(a) and — P<j e (a), decrease monotonously 
and tend to the asymptotic value (w e — Wo)p d ' e /(l + wq), 
the sign of which depends on the ratio of values of w e 
and Wq. It means that Pde (a) and pde{a) can change 
the sign at different moments of time depending on 
the relation between w e and Wq. For example, in the 
case of quintessential scalar field (w e , wq > —1) with 
wo > w e the pressure changes the sign from minus 
to plus at »(p =0 ) = [(1 + w )w e /(w e - w )] 1/3(1+u,c) 
and the density changes the sign from plus to minus at 
a (p=0 ) = [(1 + w Q )/(w e - wo)} 1/3(1+w ^- In the opposite 
case (w e > wo) their signs (P < 0, p > 0) always remain 
unaltered. 

Note that substituting §5§ into r.h.s. of ^ we see 
that the sign of w' is completely determined by the 
sign of — (1 + w )(w e — wo). So, for quintessential field 
w(a) is monotonously decreasing (w' < 0) function when 
w e > wq. The repulsion property of such field raises with 
time. When w e < wq, w(a) is monotonously increasing 
(w' > 0) function. The repulsion property of the field 
in this case recedes with time. For the phantom field, 
for which wq < — 1, the conditions for raising and reced- 
ing repulsion properties of the field are opposite. In the 
particular case w e = u>o the EoS parameter is constant, 
w' = 0. Another particular case w e = — 1 or wq = — 1 is 
simply A-term (w — const = —1). It must be noted also, 
that w(a) in the form Q does not allow the phantom 
divide (w = —1) crossing since in this point P de and pde 
become 0, but it is applicable in both ranges w > — 1 and 
w < — 1 separately. 
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TABLE I: The best fitting values and lcr confidential ranges 
of cosmological parameters in the CSF+CDM model deter- 
mined by the Markov Chain Monte Carlo technique using the 
available observational data (for details see 18]). First col- 
umn - model of DE with decreasing EoS parameter, second 
one - model of DE with increasing EoS parameter. The cur- 
rent Hubble parameter Ho is in units kms -1 Mpc -1 , the age 
of the Universe to is given in Giga years. 
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The goal of the paper is analysis of the possibility of 
distinguishing the scalar fields with receding and raising 
repulsion by available observational data. 

The discussed above three parameters of a scalar field 
- fide, w e and wq - are quite enough for the analysis of 
background dynamics, luminosity distance - redshift and 
angular diameter distance - redshift relations. But scalar 
fields are gravitationally unstable and their clustering 
properties depend also on one more of their character- 
istic - the effective sound speed c 2 s = SP de /5p de . Indeed, 
the equations for evolution of Fourier amplitudes of den- 
sity <5(rf e ) = Spde/pde and velocity V^e) perturbations of 
scalar field in the synchronous gauge are 



\de) + 3(c^ - w)aH5( de ) + (1 + w) : 



+ (l + w) 



v (de) = 0, 



(6) 



e) + aH(l - ic 2 s )V (de) - Yf^ S (de) = 0, (7) 



0.0 




0.2 0.0 



FIG. 1: Two-dimensional mean likelihood distribution in the 
plane w e - w for the combined dataset WMAP7+SDSS LRG 
DR7+SN Union2+HST+BBN. Solid lines correspond to la 
and 2cr confidence contours. 



where k is the wave number, h = h^Js the trace of scalar 
perturbations of metric (see e.g. |18l [l9| or for gauge- 
invariant approach e.g. [13, H3 and citing therein). The 
equations for the rest of components (non-relativistic and 
relativistic) are the same as in [22j . For determination 
of the fourth parameter of dark energy, c 2 s , we specify 
the scalar field Lagrangian. We consider the scalar fields 
with Klein-Gordon (KG) and Dirac-Born-Infeld (DBI) 
Lagrangians, also called classical and tachyonic respec- 
tively, 



Ldas — X — U(4>), L tac h 



-U(0Vl-2X, (8) 



where U(<fi) and £/(£) are the field potentials, X = 
: i0 ;l /2 and X = £ ; i£'V2 are kinetic terms. For the ho- 
mogeneous background the field variables are connected 
with the variables of fluid approach by simple relations: 



U(a) = ^p de ( a ), U(a) = Pde (a)^wjaj (9) 
X(a) 



l+w(a) i x y, x 1 + w{a) 

— ^Pde{a), X(a) = . (10) 



The effective (rest-frame) sound speed c 2 ,^^ for the 



i (de) 

scalar field with given Lagrangian is defined as 

2 IU = L x 

Cs(de) - p,x 2XL,xx + L tX - 



(11) 



It equals 1 for the classical scalar field and — w for the 
tachyonic one. Therefore, these fields are quite differ- 
ent, can affect the evolution of dark matter perturba- 
tions and leave "fingerprints" in the large scale structure 
of the Universe, as it was shown in our previous papers 
[IB US IHI • The goal of this paper is the analysis of pos- 
sibility of distinguishing the scalar fields with KG and 
DBI Lagrangians by available observational data. 

For calculation of the dynamics of expansion of the 
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FIG. 2: The power spectra of CMB temperature fluctuations (left) and matter density ones (right) for cosmological models 
with classical and tachyonic scalar fields (CSF and TSF respectively) with two sets of the best fitting parameters from Table 
U The corresponding observational data from WMAP7 and SDSS LRG DR7 are shown by the filled circles. 



Universe, the evolution of perturbations in all compo- 
nents, the power spectra of matter density perturbations 
and CMB anisotropy we have used the publicly available 
code CAMB [gjj [2J| , modified to include the presented 
here expressions ([3]), (gj), (fTTj) and equations @-([7|). 



B. Determination of DE and cosmological 
parameters 

Therefore, the scalar field model of dark energy, de- 
scribed in the previous subsection, involves three pa- 
rameters Vide, w e and wq which must be determined by 
comparison of the calculated predictions on dynamics 
of expansion and large scale structure of the Universe 
with corresponding observational data. Since all predic- 
tions and data are related with other components (dark 
matter, baryons, thermal cosmic radiation) it should be 
done jointly with other cosmological parameter. Con- 
centrating on the analysis of possibility of determination 
of dark energy parameters we consider the cosmological 
model with minimal set of six parameters: density pa- 
rameter of baryons density parameter of cold dark 
matter u> c dm, Hubble constant Hq, spectral index of ini- 
tial matter density power spectrum n s , amplitude of ini- 
tial matter density power spectrum A s and reionization 
optical depth T re i- So, we have nine unknown parame- 
ters, but the number of independent ones is eight, since 
we have assumed the spatial flatness of the Universe. 
Indeed, the dark energy density parameter fide in this 
case can be calculated from the zero curvature condition: 

fide = 1 — — flcdm, where fib = Wfch , £l c dm = ^cdrrM , 

where h = 7Jo/100km/s ■ Mpc. 

To determine the best fitting values and confidential 
ranges of the scalar field parameters to geth er with other 
cosmological ones in our previous work |18| we have per- 
formed the Markov Chain Monte Carlo (MCMC) analysis 
for the set of current observational data, which include 
the power spectra from WMAP7 [H [26| and SDSS DR7 



[27j . the Hubble constant measurements [28| . the light 
curves of SN la [29[ and Big Bang nucleosynthesis (BBN) 
prior [H. 

We have used the publicly available package CosmoMC 
[Hill, which includes the code CAMB [Ulli for cal- 
culation of the model predictions for sampled sets of 8 
cosmological parameters: w e , wq, Wb, uj c d, n , Hq, n s , A s 
and T r ei- The CosmoMC code has been modified to be 
run with the proposed here parametrization of dark en- 
ergy EoS parameter ([3]) . The flat priors — 1 < wq , < 
have been used to take into account the quintessential 
properties of scalar fields with KG and DBI Lagrangians 
(lower limit) and the constraints following from observa- 
tional data related to the recombination and nucleosyn- 
thesis epochs (upper one). 

We have performed two MCMC runs for the eight- 
parametric flat cosmological model with the classical 
scalar field. Each run had 8 chains and the number of 
samples in each chain was ~ 200 000. For the first run 
we have used only the mentioned above flat prior for w e - 
The set of best fitting parameters determined by this run 
is marked by and presented in the Table Q] together 
with la confidential ranges. All parameters, excluding 
w e , are well constrained, the one-dimensional posterior 
and mean likelihood distributions are close and similar 
to Gaussian (half-Gaussian for wq), confidential ranges 
are narrow. The EoS parameter at early epoch w e is es- 
sentially unconstrained: its la confidence range is wide 
and coincides practically with the allowed by prior range 
[-1, 0]. The mean likelihood and posterior are different 
and the likelihood is bimodal, as it is shown in Fig. [T] 
The first peak is close to -1, another one to 0. The best 
fitting value of w e in the set p 1 corresponds to the first 
peak. In this case w e < wq which means that the best 
fitting scalar field model of dark energy has receding re- 
pulsion properties (w/ > 0). 

In order to obtain the best fitting parameters corre- 
sponding to the second peak of the likelihood distribu- 
tion we have performed analogical run with additional 
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condition w e > wq. The set of best fitting parameters 
determined by this run is marked by p 2 and presented 
in the Table HI Now, the best fitting value of w e cor- 
responds to the second visible in the upper left corner 
of Fig. [T] peak of mean likelihood distribution. In this 
case w e > wo and we can say that the best fitting scalar 
field model of dark energy has raising repulsion proper- 
ties (w' < 0). 

As we see, large variation of w e does not change essen- 
tially other parameters: each parameter from the set p 2 
is in the la range of the corresponding one from the set 
p-L and vice versa. The — logZ's (last row of the TableHJ 
for both sets are very close. The power spectra of matter 
density perturbations and CMB temperature fluctuations 
for both sets of the best fitting parameters p x and p 2 are 
presented in Fig. [5] For the same sets of parameters we 
have calculated also the spectra for tachyonic scalar field 
models of dark energy. For all four models corresponding 
lines are superimposed and well match the observational 
spectra. Therefore, we have double degeneracy: in type 
of the dynamics of scalar field (receding-raising repulsion 
properties) and its Lagrangian (classical-tachyonic) . Is it 
possible to distinguish them in principle? It is subject of 
the next sections. 
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FIG. 3: The dependences of EoS parameter (top panel) and 
ratio pde/pm (bottom panel) on the scale factor for models 
with best fitting parameters pj and p 2 . The current epoch 
corresponds to a = 1. 



III. DYNAMICS OF EXPANSION OF THE 
UNIVERSE, SN IA LIGHT CURVES AND 
STATEFINDERS 

At first we analyze the evolution of scalar field param- 
eters in the homogeneous (background) Universe. The 
time dependences of the EoS parameter and the ratio 
of scalar field energy density to dark matter density in 
the models with best fitting parameters p ± and p 2 are 
presented in Fig. [3] As we see, in spite of closeness 
of these values for both parameter sets in the current 
epoch (a = 1), they are quite different in past (a < 1) 
and future (a > I). The EoS parameter in the model 
with p-L increases monotonically from -0.97 at the be- 
ginning of expansion to -0.93 at the current epoch and 
will continue increasing in future up to discontinuity of 
the second kind, caused by the energy density crossing of 
zero (a(p dc =o) = 950). Immediately after that the Uni- 
verse reaches the turnaround point (a ta = a( Pde=0 ) + Sa 7 
8a /a <C 1) and begins to collapse 18]. At the early 
epoch the scalar field energy density is insignificant simi- 
larly to the cosmological constant. In the model with p 2 
the EoS parameter decreases monotonically from -0.05 
at the early epoch to -0.99 now and will asymptotically 
approach -1 in far future. The energy density in the 
early epoch traces the density of dark matter. In fu- 
ture it will asymptotically approach the constant value 

Pde^ — Pd°e( w e ~ u 'o)/(l + w e ). So, the future of the 
Universe with such field is de Sitter expansion. 

The evolution of potentials and kinetic terms of clas- 
sical and tachyonic scalar fields for models with the best 
fitting parameters p 1 and p 2 is shown in the Fig. 0J One 



can see, that both fields are quite different for p x and 
p 2 parameter sets, but for each of them the potentials of 
both fields are similar, while kinetic terms are different. 
In all cases scalar fields roll slowly to the minima of their 
potentials. In the models with p : parameter set it will be 
reached in finite time (« 268 Gyrs, turnaround point), 
in the models with p 2 one - at infinite time. Near the 
turnaround point in the model with p 1 the kinetic terms 
of both fields become dominating (top panels of Fig. 0} , 
quintessential scalar fields behave as the k-essential ones. 

Using the dependences of densities of each compo- 
nent on the scale factor one can deduce from the Einstein 
equations the following equations for the background dy- 
namics: 

H = lWa> 4 + ^™/« 3 + JW(a), (12) 
= 1 2n r /a 4 + n m /a 3 + (1 + 3w)n de f(a) 
2 flr/a 4 + fl m /a 3 + ttdef(a) 

where f{a) = [(1 + w )a-^ 1+w ^ + w e - w ]/(l + w e ), 
Q m — fli, + fl cc i m is the density parameter of matter (non- 
relativistic) component at current epoch, Sl r = f2 7 + 0„ 
is the density parameter of relativistic component at cur- 
rent epoch, H = a/a is the Hubble parameter (expansion 
rate) for any moment of time and q = —ii/aH 2 is the 
acceleration parameter for any moment of time. They 
describe completely the dynamics of expansion of the ho- 
mogeneous isotropic Universe. 

The temporal dependences of H and q for models with 
p : and p 2 are shown in Fig. [SJ It can be seen that their 
past evolution (a < 1) is practically indistinguishable, 
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FIG. 4: Evolution of potentials and kinetic terms of classical (left column) and tachyonic (right column) scalar fields for models 
with best fitting parameters p 1 (top panel) and p 2 (bottom panel). 



but future one is different. In both models the accelerated 
expansion began at a « 0.58. At current epoch the accel- 
eration parameters are qo ~ —0.57 in the model with p 1 
and go w —0.52 in the model with p 2 . In the last model it 
will continue decreasing in future, approaching asymptot- 
ically -1 as in de Sitter inflation. In the model with p x the 
acceleration parameter will reach minimal value q m in ~ 
—0.87 at ct(g min ) ~ 4.95 and then will begin to increase. 

At a {q=0) = [(1 + w )(l + 3w e )/(w e - Wo )/2] 1/3(1+ ^ a 
570 the accelerated expansion will be altered by de- 
celerated one. The pressure will become positive at 

i 

O(u;=o) > [w e (l + wo)/(w e — wo)] 3 (!+™e) a 678 and in 
finite time according to (|5|) will reach the constant pos- 
itive value fj™ = ( w o — w e)p^e ~ 0-035p^ causing 
the turn around at a » 960 and collapse. So, the cosmo- 
logical model with p x has future finite-time Big Crunch 
singularity. 

Fig. [5] illustrates the model degeneracy: two different 
scalar field models of dark energy are indistinguishable 
by the dynamics of expansion of the Universe. How deep 
is this degeneracy? To answer this question we have cal- 
culated the relative differences AH(a)/H(a) and Aq/q 
of models with p 1 and p 2 . The results are presented in 
Fig. |51 One can see that the appreciable difference (1%) 
of H in models with p : and p 2 in the past (a < 1) ap- 
pears at redshifts z > 1 (a < 0.5) and reaches 2% at 
high z. The relative difference of acceleration parame- 
ters in these models is maximal (« 8%) at current epoch 
(a « 1) and rs 5% at the beginning of accelerated expan- 
sion (a » 0.4 - 0.6). 

Since both values are usually determined from the lu- 
minosity distance-redshift relation, it is interesting to 



compare the difference of distance moduli 

(m-M) = 5 log d L + 25, (14) 

where 

c , , f z dz' 



d L = 77-(l + 2) 

H ° Jo ^Jn m (i + z'Y + n de f(j^r) 

for these models with observational uncertainties of their 
determinations from the light curves of SN. In Fig. [7] we 
present the distance moduli for both models and rela- 
tive difference of them together with corresponding cur- 
rent observational data [3l|. Both models match obser- 
vational data on SN SDSS distance moduli equally good, 
the difference between them does not exceed 0.1%, which 
is essentially lower than dispersion of observational points 
around of the best fitting curve. To distinguish them 
other tests based on accurate measurements of accelera- 
tion parameter in the vicinity of the Local Group should 
be proposed or radical improvement of the accuracy of 
existing ones should be made. 

If in the future high-precision observations it is 
possible to measure the second derivative of H then 
statefinder parameters r = 1 + 3H/H 2 + H/H 3 and 
s = (r — l)/[3(<7 — 1/2)], introduced by (34J, can be used 
for distinguishing the scalar field models with different 
time variable EoS parameters. For matter plus dark 
energy dominated epoch they can be presented in our 
paramctrization ([3]) as follows: 

r = 1 + 4.5(1 + w)w e flde(a), s = (1 + w)w e /w, 

where flde{<i) = 8wGpde(a)/5H 2 . Their dependences on 
scale factor a for the best fitting models p x and p 2 are 
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FIG. 5: The dynamics of expansion of the Universe with 
barotropic dark energy: H(a) (top panel) and q(a) (bottom 
one) for the pj and p 2 models. 
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FIG. 6: The relative differences AH(a)/H(a) and Aq/go of 
models with p 1 and p 2 . 



FIG. 7: Top panel: The distance moduli (m - M)(z) for 
models with p 1 and p 2 (lines) and SN SDSS data with 
different fittings of SN la light curves (signs). Bottom 
panel: the relative difference A(m — M)/(m — M) of two 
models (line) and scattering of data points from the top 
panel around the model line with best fitting parameters p^ 
[(m - M)obs - (m - M) mod ] /{m- M) mod . 
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shown in Fig. [8] One can see that the differences between 
s(p 1 ) and s(p 2 ) at high z as well as between r(p 1 ) and 
r(p 2 ) at low z are essentially larger than for parameters 
H and q (Fig. [5] and©. 

IV. MATTER DENSITY AND CMB POWER 
SPECTRA 

Let us analyze the large scale structure and CMB 
anisotropy characteristics calculated for models with 
p 1 and p 2 parameter sets, their differences and model 
distinguishability in the light of them. The dark energy 



FIG. 8: Statefinder parameters r and s for models with p : 
and p 2 . 



affect these characteristics via dynamics of expansion of 
the homogeneous Universe (scale-independent growth 
factor) and gravitational interaction between pertur- 
bations of all components for the same fc-mode. The 
evolution of scalar field density perturbations depends 
also on Cg, therefore the influence of them on the 
matter density ones will be different for KG and DBI 
Lagrangians. Using the modified CAMB code we have 
integrated the Einstein-Boltzmann system of equations 
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for all components (scalar field, dark matter, baryons, 
thermal radiation, massless active neutrinos) in the 
cosmological models with parameter sets p x and p 2 . 
The results for baryons, cold dark matter and scalar 
fields with KG/DBI Lagrangians (wave number of 
perturbations k — 0.05 Mpc -1 ) are shown in Fig. [9] 
In the case of model with p x the evolution tracks of 
scalar field density perturbations are similar for both 
Lagrangians, but they differ in the case of model with 
p 2 . Since scalar field density perturbations are not 
observable, let us look how they influence the evolution 
tracks of matter density perturbations. For this purpose 
in Fig. [10] we present the differences £±5( c dm)l &{cdm) — 
(£( c dm)(Pi;a) - S {cdm) (p 2 ;a)) /<5( c(im )(p 1 ; a) 
for KG and DBI scalar field models 
(left panel) as well as A5( cdm )/6( cdm ) = 
{S(cd m ) (CSF; a) - (5 (crJm )(TSF; a)) /S {cdm )(CSF; a) for p x 
and p 2 parameter sets (right panel). The differences 
between models with p 1 and p 2 for scalar fields with KG 
and DBI Lagrangians reach ~ 2% at the current epoch. 
In the case of p x parameter set the scalar field with 
DBI Lagrangian is indistinguishable from the scalar field 
with KG one (dashed line in right panel of Fig. flTJ|) . In 
the case of p 2 parameter set the difference between both 
fields does not exceed 0.3%. It means that at the current 
level of accuracy of observational data and numerical 
codes for calculations of model predictions they cannot 
be distinguished too. For final conclusion other fc-modes 
should be analyzed. 

For such purpose we have calculated by the modified 
version of CAMB the power spectra of matter density 
perturbations for the same four models P(CSF, p x ; k), 
P(CSF,p 2 ;k), P(TSF, Pl ;k) and P(TSF,p 2 ;k) (right 
panel of Fig. [5]). The relative differences AP/P are 
presented for them in Fig. [TTJ The observational relative 
la errors of SDSS LRG DR7 data [13| are shown there 
too. The maximal differences between spectra in models 
with p x and p 2 are at k ~ 0.01 h/Mpc (w 8% in the 
models with CSF and « 6% in the models with TSF). At 
lower scales, k ~ 0.1 h/Mpc, they are ~ 4 — 5%, that is 
lower than the observational uncertainties at this scales. 
In the case of p x parameter set the models with CSF 
and TSF are indistinguishable, the relative difference 
is less than 0.1%. In the case of p 2 parameter set it is 
maximal (» 1.5%) at k -0.01 h/Mpc, but it is « 0.6% 
at k —0.1 h/Mpc, where observational uncertainties 
are lowest, « 6%. So, the possibility of distinguishing 
CSF from TSF by P is unpredictable. And vice versa, 
distinguishing of the model with p 1 from one with p 2 
by the matter density power spectrum data expected in 
upcoming decade looks possible. 

Let us analyze now the possibility of distinguishing 
these four models by observational data on CMB 
anisotropy. In Fig. [2] the angular power spectra of 
CMB temperature fluctuations Cj T [12[ for them are 
presented in left panel. All four lines are visually 
superimposed and go near most points in their la error 
bars. The relative differences between them ACj T /Cj T 



as well as relative la errors of WMAP7 data [25|, [26| are 
shown in Fig. [TJ] The maximal differences (~ 3 — 4%) 
between Cj T of models with p 1 and p 2 are at low 
spherical harmonics, where cosmic variance is too large 
to distinguish between such models. In the range of 
acoustic peaks £ — 100 — 700, where observational data 
are most accurate (— 2%), the difference between Cj T 
of models with p 1 and p 2 is somewhat smaller (< 1.5%). 
So, the expected data releases on CMB anisotropy 
from WMAP and Planck teams would probably allow 
the possibility to answer the question "which scalar 
field dark energy, with decreasing or increasing EoS 
parameter, fills our Universe?" . To answer the question 
"which is the field Lagrangian?" it will be harder in 
the case of decreasing EoS parameter and practically 
impossible in the case of increasing one (see right panels 
of Figs. [Mill). 

To support the latter conclusions let us compare 
the relative differences of CMB temperature fluctua- 
tions ACf T /Cf T and polarization ACf E /Cf E power 
spectra in models with both fields and parameter sets 
with the observational uncertainties, modelled for the 
Planck satellite in the following way. Assuming that 
the noise part is due to the combined effect of Gaussian 
beam and spatially uniform Gaussian white noise, for 
the experiment with known beam width and sensitivity 
the noise power spectrum for each channel can be 
approximated as follows: 



= 9 2 fwhm rf exp 



£{l+l) 



9 2 

fwhm 

8 log 2 



where j stands for either TT or EE, 9t W h m is the full 
width at half maximum of the Gaussian beam and cr, 
is the root mean square of the instrumental noise. The 
non-diagonal noise terms vanish since the noise contri- 
butions from different maps do not correlate. For the 
experiments with more than one channel the total noise 
power spectrum is obtained as: 



1 n ehan 

i = V 

AT jj(tot) A^/ 
iv / »=1 



1 



N, 



n(i) ' 



where n c han is the number of channels. The described 
procedure was proposed by [38[ and implemented in their 
code FuturCMB [39|], which we use here. 

In Fig. [13] we show the estimated errors for the Planck 
experiment with 3 channels (for each of them 0f w hmi a T 
and aE are 9.5 arcmin, 6.8 fiK per pixel and 10.9 \xK 
per pixel; 7.1 arcmin, 6.0 fiK per pixel and 11.4 fiK per 
pixel; 5.0 arcmin, 13.1 /iK per pixel and 26.7 fiK per 
pixel correspondingly). The observed sky fraction is as- 
sumed to be f s ky = 0.65. The models with the same 
fields but different parameter sets p x and p 2 can be dis- 
tinguished by the data with such precision. For the CMB 
temperature fluctuations power spectrum the difference 
between studied models exceeds the estimated error level 
at high spherical harmonics, while for the polarization 
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FIG. 9: Evolution of linear density perturbations of cold dark matter (dash-dotted line), baryons (dotted) and scalar field 
(solid) for models with p 1 and p 2 . Classical scalar field - left column, tachyonic one - right. The wave number of perturbations 
is k = 0.05 Mpc" 1 . 




FIG. 10: The influence of scalar field perturbations on the evolution of linear matter density ones (k = 0.05 
Mpc -1 ). Left panel: the relative difference of CDM density perturbations in models with sets of parameters p x and p 2 
A(5( cdm )/5( c dm) = (<5(cdm)(Pi;a) - <5( cdm )(p 2 ;a))/5( cdm )(p 1 ;a) for CSF (solid line) and TSF (dashed line). Right panel: the 
relative difference of CDM density perturbations in the models with CSF and TSF A5( cdm )/(5( ct jm) = (5( c dm) (TSF; a) — 
<5( c dm)(CSF; a))/5( ct j m )(CSF; a) for the sets of parameters p 1 (solid line) and p 2 (dashed line). 



power spectrum at low spherical harmonics, where it is 
maximal. The models with different fields but the same 
parameter sets (corresponding to both decreasing and in- 
creasing EoS parameters) are still indistinguishable at 
such level of experimental precision. 



V. DISCUSSION AND CONCLUSION 

Weakness of constraints on dark energy parameters 
for models with time variable EoS is widely discussed 
in the literature [35H37j. Moreover, the degeneracy be- 
tween them limits strongly the possibility to test whether 
w is constant or not. Our results, presented in the Tab. 
2 of [llj], confirm that. We have used the set of obser- 
vational data including the light curves of SN la, the 
power spectra from WMAP7 and SDSS DR7, the Hub- 
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FIG. 11: Left panel: the relative difference of matter density power spectra AP/P in the models with best fitting parameters 
Pj and p 2 for classical and tachyonic scalar fields. Right panel: the relative difference of matter density power spectra AP/P 
in the models with classical and tachyonic scalar fields for two sets of the best fitting parameters p 1 and p 2 . Dots show 
observational uncertainties (la) of SDSS LRG DR7 data. 




FIG. 12: Left panel: the relative difference of CMB temperature fluctuations power spectra ACj T /Cj T in the models with best 
fitting parameters p 1 and p 2 for classical and tachyonic scalar fields. Right panel: the relative difference of CMB temperature 
fluctuations power spectra ACj T /Cj T in the models with classical and tachyonic scalar fields for two sets of the best fitting 
parameters p 1 and p 2 . Dots show observational uncertainties (la) of WMAP7 data. 



ble constant measurements and BBN prior to analyze the 
possibility of determining the parameters of dynamical 
scalar field jointly with other cosmological ones. Pro- 
cessing of the MCMC chains for eight-parametric mod- 
els (p = (fi de , w , w e , oj b , uj cdm , H Q , n s , A s , T re i)) gives 
marginalized posterior and likelihood distributions as 
well as the best fitting values of parameters and their 
confidential ranges. For the most parameters the poste- 
rior and likelihood distributions are unimodal and sim- 
ilar, the confidential ranges are narrow. The exception 
is w e , for which the marginalized posterior distribution 
is unimodal, but marginalized mean likelihood distribu- 
tion is bimodal. The values of maximal likelihood for 
both peaks are very close. The first peak corresponds 
to the scalar field model of dark energy with increasing 
EoS parameter, the second to model with decreasing one. 



The question is: which accuracy of observational data is 
necessary to distinguish these two models of dark energy? 
We have shown, that the difference of luminosity distance 
moduli A(m — M) in the models with decreasing and in- 
creasing EoS parameter does not exceed 0.1% (Fig. [7]), 
while the most accurate observational data on SN lumi- 
nosity distances (SN SDSS, MLCS2k2 and SALT2 [U) 
disperse around the best fitting model line in the range of 
~ 2%. Therefore, the radical improvement of statistical 
and systematic uncertainties of data is necessary for dis- 
tinguishing between these models of dark energy. Then 
the statefinder parameters r and s 3J] allow the possi- 
bility to differentiate effectively the scalar field models of 
dark energy with increasing and decreasing EoS param- 
eters. 

To diminish or remove the degeneracy between dark 
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FIG. 13: Left: the relative differences of CMB temperature fluctuations ACj T /Cj T and polarization AC EE /C EE power 
spectra (from top to bottom) in the models with classical and tachyonic scalar fields for two sets of the best fitting parameters 
p x and p 2 . Right: the relative differences of CMB temperature fluctuations ACf T /Cf T and polarization AC EE /Cf E power 
spectra (from top to bottom) in the models with two sets of the best fitting parameters p 1 and p 2 for either classical or 
tachyonic scalar field. Dots show the observational uncertainties, modelled for the Planck satellite. 



energy, dark matter and curvature parameters other data 
must be also improved. The difference of the mat- 
ter power spectra is maximal (~6-8%) at scales 0.1-0.5 
/i _1 Gpc (left panel of Fig. ITT]) . The maximal accu- 
racy of matter density perturbations power spectra ob- 
tained in the current sky surveys reaches ~6% at scales 
~ 50 /i _1 Mpc (SDSS DR7), which is somewhat above the 
difference between model power spectra at this scales. If 
expected data from the next releases are 1.5-2 times more 
accurate, then they allow us to distinguish surely the 
models with decreasing and increasing EoS parameters. 
The important data in the used set are the data on CMB 
anisotropy, as the difference of angular power spectra of 
CMB temperature fluctuations is ~l-4%, which is some- 
what smaller than accuracy of the current measurements 
by WMAP at scales - 0°.5-2° (i ~ 100-500) (left panel 
of Fig. [T2|) . We believe that the expected data release 
from WMAP team (nine-year observations) and particu- 
larly expected Planck data together with other improved 
data will make it possible to determine whether the EoS 
parameter of dark energy is increasing, decreasing or con- 
stant. 

Possibility of distinguishing the scalar field models of 
dark energy with different Lagrangians is less optimistic. 
We have analyzed only two from many proposed in the 
literature ones and found that differences between pre- 
dictions of such models are essentially lower than accu- 



racy of current observational data. So, the scalar fields 
with KG and DBI Lagrangians are indistinguishable at 
the current level of accuracy of observational data (right 
panels of Fig. [T0TIT2^) . Only in the case of decreas- 
ing EoS parameter the scalar fields with KG and DBI 
Lagrangians are potentially distinguishable, perhaps al- 
ready by Planck data. But in the case of increasing EoS 
parameter they are practically indistinguishable, since 
differences for them are <C 0.1%, which is essentially 
lower than the accuracy of theoretical and computer cal- 
culations of predictions. These small differences have the 
simple explanation: in such case the densities of scalar 
fields were too low in the past and did not leave the 
fingerprints in the CMB anisotropy and LSS. This con- 
clusion can be generalized for any type of scalar fields - 
quintessence, phantom, quintom, etc. 

We have considered the principal possibility of distin- 
guishing the scalar fields with decreasing and increas- 
ing EoS parameters and two Lagrangians in cosmological 
models with the minimal set of parameters. Adding extra 
parameters - masses of active neutrinos, sterile neutrino, 
tensor mode, primordial magnetic fields, etc. - weakens 
constraints on dark energy parameters and possibility to 
distinguish the different scalar fields modelling it. So, 
only join progress of theory, laboratory experiments, as- 
trophysical and cosmological observations will lead to the 
unveiling of nature of dark energy. 
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